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Abstract 

In this paper we introduce a simple phenomenological model of the conduction be- 
tween a couple of serial or parallel quantum dots. This model is extended to an arbitrary 
of number and to a square array of quantum dots. The local potential is computed taking 
'-^ into account the net charge at the quantum dot. Master equations are presented for total 

current and charge number N considering both the local potential and the local current. 
Interesting results are reported, namely the Negative Differential Resistance for serial 
quantum dots as well as resonant conductance. In parallel configuration two independent 
conduction channels have been observed. 
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^ 1 Introduction 

Single-electron devices are currently conceived to take advantage of the tunnel current between quan- 
tum states belonging to nano-scale particles [I] [2] . Up to now, the study of such tunnel current has 
been carried out using simple phenomenological (i.e. rate equation) approaches OH] or more complex 
quantum theories including NEGFF El [7] , transfer Hamiltonian approach [H] |9] . Concerning the 
first case, simple rate equations give us useful and intuitive approaches. The effects of local potential 
due to selfcharge can be also included by means of a couple of master equations [3]. On the other hand, 
quantum study of non-equilibrium schemes results in a more complex approach. 
Until now research has been mostly concentrated on single dots, but rapid progress in microfabrication 
technology has made possible the extension to couple-dot system with aligned levels [10], [11] [12] . A 
very interesting result was obtained by Shanguan et al |13j . These authors describe the transport in 
one-dimensional quantum dot array by using many body approaches. Up to now the only computation 
of transport in an array of quantum dots was done by Carreras et al [14] but no local potential due 
to selfcharge was included. 

In this work we present a semiclassical approach based on non-coherent rate equations [U [15] that 
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is extended from two to an arbitrary number of quantum dots either serial, parallel and square ar- 
ray. The whole array of quantum dots has been treated as a separated quantum system. Then the 
leads (reservoirs) were incorporated through the rate equations for single dots. In this approach, the 
coherent effects does not appear and the probabilities of finding an electron in two different dots are 
coupled by some rate 7. This treatment is called classical (non-coherent) approach. 
The effects of the local potential are computed within the selfconsistent field (SCF) regime. The SCF 
method is widely used because the exact method based on a multielectron picture is usually impossible 
to implement. As example we compute the total current and charge transport distribution iVj for two 
parallel and serial quantum dots (figure 1). The model has been implemented in a home made pro- 
gram in MATLAB code that gives us the current voltage characteristics, the differential conductance 
and the population inside each dot for a range of applied voltage. Some known results as Negative 
Differential Resistance (NDR) and resonant conductance are well explained. Moreover we found that 
the parallel quantum dots act like two independent conduction channels. 




Figure 1: System composed of quantum dot 1 and 2, and two leads. % are the tunneling rates 
between each part of the system, /i^ and fi^ are the electrochemical potentials in right and 
left lead, respectively. Representation of the two different systems: serial configuration (a) and 
parallel configuration (b). 



2 Theroretical model 

When we have two leads coupled to a system with no voltage applied between them, the system is in 
equilibrium with a common electrochemical potential. In this equilibrium state, the average number 
of electrons in any energy level is given by the Fermi function: 

f{E) = l + exp((E- fi)/k B T) (1) 

An applied bias voltage Vq = fiL ~ >MR drives the system out of equilibrium. Each contact tries to bring 
the system into equilibrium with itself. The system is thus forced into a balancing act between two 
reservoirs with different population. Then, the average number of electrons N inside the system will 
be something between fn and fi. Rate-equation calculation consider the dynamic behavior of the 
number of electron rii at each quantum dot balanced with the incoming and outcoming flux either 
from a partially filled (/l,_r) lead or other quantum dot (rij) [TBI IT7]. 
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2.1 Serial quantum dot 

Rate-equation for serial Qd can be written as follow: 

-jjT = JUL + 72«2 - (71 + 12)11! (2) 

^ = 73/r + 72«i - (72 + 73)"2 (3) 

where coefficients ji are the transition probabilities as described in figure la. Sl and Sr are the Fermi- 
Dirac distribution functions in the electrodes. Density of states at electrodes is considered constant 
in the energy range that transport takes place, rij is the average number of electrons at i — th Qd. 
The steady state solution of (2) and (3) give as the expression 

71 (72 + 73) h + 7273/R u \ 

U\ = — 7y (4 J 

(71 + 72) (72 + 73) - 72 

73(72 + 7i)/fl + 7271 h 
(71 + 72) (72 + 73) - 72 
Besides the out coming net flux from Qd2 to the right lead is 



6\ 2 T I LJ J H 1 12 LJ L 

n 2 = 77 , \ , ~ \ 72 ( 5 ) 



I = (-q)^(fR-n 2 ) (6) 
That help us to consider a n-independent way to express the current using (4) and (5) 

I s = Ul ~ Sr) (7) 

7i 72 + 71 73 + 7273 

The expression of the current can be physically interpreted quite easily: the three barriers in the Qd 
stack act as a series conection of resistors corresponding to the inverse of tunnelling rates Rtotai = 

2.2 Parallel quantum dot 

We can propose the same procedure for parallel Qd. In this case rate-equations are 

^jjr = 7i/l + 73/i? + 72"-2 - (71 + 72 + 73)«i (8) 

^jjr = 7i/l + 73/i? + 72m - (71 + 72 + 73)^2 (9) 
considering the figure lb. The steady-state solution of (8) and (9) is 

7i(7i + 73 + 27 2 )/l + 73(71 + 73 + 272) Sr /irA 
n\ = -. r-^ — I1UI 

(7i + 73) + 27273 + 27271 

n = 7i(7i + 73 + 27 2 )/l + 73(71 + 73 + 272)//; , . 

(71 + 73) 2 + 27273 + 27271 

The net out coming flow is 

I P = (-Q)^(SR-ni) (12) 

for each quantum dot. And substituting (10) and (11) in (12) we obtain the n-independent expression 
of the current 

T p 7173(71 +73 + 272) ,. , v n „x 

1 7 1 \2 1 o Ul-jr) (13) 

(71 + 73) + 27372 + 27271 
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2.3 Effect of broadening 



Up to now we are considering only one energy level per dot. But we have missed the broadening of 
the level that inevitably accompanies any process of coupling to it. The standard way to introduce 
energy level broadening as a consequence of contacts is to assign a Lorentzian shape to density of 
states [TBI H2| 

The broadening 7 is proportional to the strength of the coupling. 7 = J2i li , where 7, are the 
transition rates introduced before. This comes out of a full quantum mechanical treatment, but we 
could rationalize it as a consequence of the "uncertainty principle" that requires the product of the 
lifetime (= ft/7) of a state and its spread in energy (7) to equal Ti. The current expression is modified 
into: 

I = (-q)^(fR ~ ni) -+ I = (-<?) J j^PRPi(E)(f R - m)dE (15) 

where now m can be interpreted as the non-equilibrium distribution functions. This allow us to 
compute the total number of electrons at each Qd as 



/+00 
dEpim i = l,2 (16) 
-00 



This procedure allow to rewrite n-independent equations as 

/' = ! f 00 dE 717273/WLP1P2 ^ _ ( 

jP = Q r g 3 H7z{llpL + l3pR + l2{Pl+ P2))pLpRpi „ ^ 
1 hJ- 00 (j 1 p L +J 3 p R ) 2 +y 3 j2PR(Pl+P2)-hl72PL{Pl+P2) 

for serial and parallel currents. Where pn and pi are the density of states in the leads (constants), 
while p\ and P2 are the DOS in each Qd. 



3 Local potential 



Another factor that we have to take into account is how the voltage applied to the external electrodes 
change the electrostatic potential inside each dot. It is easy to see that this can play an important 
role to determining the shape of the current- volt age characteristics. The classical right solution of 
potential at each quantum dot involves Poisson equation: 

V(e r W) = -^ (19) 

Where solution is 

Q 2 

U = U L + ~ — AiV (20) 

^tot 

Being Ul the Laplace solution of the system. Considering the capacitive coupling and boundary 
conditions 

^ = E^(-^) ( 21 ) 

■ Ctot 

Where i covers all the elements of the system, including electrodes and arbitrary number of Qd, and 
Ci is the capacitive coupling between i and i + 1 element. Finally Ctot = J2i C«- The charge energy 
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system constant U$ = q 2 /Ctot gives us the potential increase as a consequence of the electron addition. 
We observe that potential depends on the increasing charge density, but at the same time charge 
density depends on energy level, a function of local potential. This considerations confers equation 
for N (16) and for U (23), the position of master equations, as usual in similar cases (5) and (6). 
For serial Qd we can consider that first Qd is capacitive coupled with left lead and second Qd, while 
the second Qd is only coupled with the first Qd and the right lead. Thus the Laplace potential results 
are 

U S L1 = ^(-qV G ) + ^U 2 (22) 

<^tot <^tot 

U S L2 = £-U x (23) 

<^tot 

where we consider Vq the potential of left electrode and for the right one. C s and C c are coupling 
capacities between Qd and electrode and between Qds. Ui is the total potential in the i — th Qd. 
Similar expressions can be derived for the parallel system: 

U[ 1 = ^(-qV G ) + ^U 2 (24) 

<^tot ^tot 

U[ 2 = ^(-qV G ) + ^ (25) 
<^tot L'tot 

We shall solve U\ and U 2 equations simultaneously in an iterative process forcing comparison of results 
[20] and thus ensuring simultaneous convergence. 



4 Results 

For simplicity we consider fully symmetric system 7$ = jj and t% = €j . The effects of local potential on 
DOS should be computed as p(E) — > p{E — U), which will modify computations for population N and 
currents /. Transport requires energy levels of Qds to be between electrochemical potentials of leads. 
Moreover overlapping of the DOS is also required. Both conditions ensure electron transport and net 
current. In this model we consider 73 << Uq and KbT « Uq. Both conditions are necessary for 
Coulomb Blockade [21]. The self-consistent charging model based on the Poisson equation represents 
a good zero-order approximation (called Hartree approximation) to the problem of electron-electron 
interactions, but it is generally recognized that it tends to overestimate the effect. Corrections for 
the so-called exchange and correlation effects are often added, but the description is still within the 
one-electron picture which assumes that a typical electron feels some average potential, U, due to 
the other electrons. Despite our self-consistent potential model do not explain potential differences 
between spin up and down as consequence of broken spin degeneration |22l 123] . we can obtain a first 
approach to Coulomb Blockade. 

Figures 2a and 2b show obtained intensities for parallel model in a range of capacities values. The 
/(V) curve is strongly dependent on C s and C c through the Laplace solution of the system. The 
value of the capacity indicates how the system is coupled. When C s tends to cero, the two Qd are 
electrically decoupled with the leads and the potential has not Laplace term, this case implies a lightly 
coupled system. The other case is when the value of the capacity between Qd (C c ) goes to cero, which 
means no electrical influence between Qd. The Poisson term follows the charge Ni at each Qd and 
it is therefore always positive, inducing a shift on the potential. Conductance is depicted in figures 
3a and 3b for the two ranges of capacities. In coherence with results shown in |22] conduction gap is 
found. The peak of conduction becomes finite as consequence of the temperature. Furthermore the 
charge effects tend to shift and widening the peak. 
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Voltage (V) 




Figure 2: I(V) characteristics calculated in the parallel configuration in a different range of 
parameters C G and Cs- C s2 is a constant . In the upper figure (a) a rectifying effect is obtained 
due to the weak coupling between the Qds and the leads. In (b) the effect of the electrostatic 
coupling between the two Qds is shown. fi R = 0,7i = 72 = 73 = 0.005e\^, K B T = 0.025el/, 
U = 0.25eV, ei = 0.2eV, e 2 = 0.2eV. 



Concerning the symmetry let us take into account two considerations. First we compute the Laplace 
domain at each Qd, as follow 

ei=e +-^(-qV G ) + ^U 2 where U 2 = ^-(-qV G ) + -^U 1 (26) 

Wot Wot Wot Wot 

that define a recurrence, (24) and (25). In order to analyze the results we will extend the series for 
the energy level of one dot. Thus energy level becomes 

ei = eo + £-(rqV G ) + ~~{-qV G ) + S%U X * e + ^^(-gF G ) (27) 

Wot Wot Wot ^tot + W2 

Where we have used that the recurrence is a geometric series, a + ar + ar 2 + . . . = The value C s2 is 
the coupling capacity with the second lead. As previously stated, the necessary condition for transport 
is energy level must stay between the electrochemical potentials \i\ > ej > fi 2 of the leads. In parallel 
Qd system, there are multiple pathways possible for electron transport and thus the condition e\ = e 2 
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Figure 3: Conductance characteristics obtained in the parallel case, in the same cases as before, 
(a) Resonant conductance peak appears as a result of the rectifying effect. When Cs increases 
the peak diminishes and the rectifying effect disappears, (b) Conductance gap is obtained . 



is not required. Then, transport starts when lead electrochemical potential reach the minimum energy 
level 6j. The potential threshold becomes: 



Cs+C. 
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Vg= g- (29) 

C s +C 3 2 

Two potential thresholds appear, for positive and negative voltage, as a consequence of previous con- 
dition. This creates a region in which no transport occurs because the previous condition is not 
accomplished, hence a conductance gap appears. The symmetry of the I(V) curves is reached when 
Vq =\Vq\. This condition occurs when eo = or C s = C S 2- In the case of a non-degenerate semicon- 
ductor the only possibility of broken symmetry is the different capacitive coupling between the two 
leads. As we can see from the expression (27) and its equivalent for 62, if the two Qds have the same €0 
and equal capacity values each Qd acts like a independent conduction channel due to the electrostatic 
potential is equal in each dot. Furthermore the electrostatic coupling between neighbours is a second 
order effect. 
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Figure 4: (a) /(V) curve of the double dot in serial configuration, showing rectification effect 
when it is weakly coupled with the lead and the two levels are aligned. The rectifying effect 
disappears and resonant peaks appears when the Qd is strongly coupled with the lead, (b) 
Resonant peaks obtained due to the different electrostatic coupling in each dot. When the 
voltage increases, the separation between the energy levels increases and the overlapping be- 
tween the DOS decreases, as a consequence, a NDR appears, hr = 0,ji = 72 = 73 = 0.005ey, 
K B T = 0.025eV, U = 0.25eV, e x = 0.2eV, e 2 = 0.2eV. 



Figures 4a and 4b show obtained intensities for serial model in a range of capacity values, in the same 
cases as before. Now, the electron has single transport way and conditions for transport are more 
restrictive. Overlapping of the density of states between Qds is also required. Conductance is depicted 
in figure 5a and 5b for the two ranges of capacities. 

As we did before, the evolution of the energy levels inside each Qd taking into account (22) and (23) 
are 

ei = + £-{-qV Q ) + £-U 2 « eo + £±£-(-qV a ) (30) 
e 2 = 60 + -^(t/i) « 60 + °" (-qV G ) (31) 

^tot + 

Where we have used the same geometric recurrence as before. For simplification we have considered 
C s = C S 2 ■ Now the two conditions imposed before in order to transport occurs and can be summarized 
as \i\ > El = E2 > /i2- As we can see from the expressions (30) and (31), ei increases faster than 
62 and thus the condition ei = 62 can not be accomplished. Introduced DOS broadening relaxes 
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Figure 5: Conductance characteristics obtained in the serial case, in the same cases as before. 
NDR and rectification effect (a). Evolution of the resonant peak, maximum of conductance, 
and the NDR as a function of the coupling capacity between dots (b). 

restrictive conditions for transport as a result of energy channel overlapping. When the voltage 
increases, the separation between ei and 62 increases as well and then the overlapping between DOS 
decrease closing some channels and decreasing the outcoming flow. As a consequence a NDR appears 
due to the different electrostatic coupling in each dot pU [25j EH EZ]- The /(V") peak is related with 
the maximum overlapping between the DOS of the first and the second Qd [11 1. Moreover the width 
is related with the DOS broadening. 



5 Conclusions 

A simple phenomenological model of the conduction between serial or parallel quantum dots has been 
developed based in rate equations for non coherent Qd system. This approach provides the most simple 
and transparent way for a description of electron transport. In the beginning we consider a system of 
two Qd linked by a ballistic channel and connected with the emitter and collector reservoirs. Level 
broadening has been introduced in the DOS changing the expression for current and the population 
inside each dot. Moreover this model can be extended easily to an array of NxM Qd. 
An important point has been the calculation of the local potential, solving the Poisson equation with 
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appropriate boundary conditions for each dot and configuration. Despite its simplicity the effect of 
self-charge has been token into account. As we may expect the calculation of the local potential inside 
each dot is the most important point. This local potential dominates de current characteristics since 
overlapping of the DOS is imposed in order to have electron transport. 

I(V) and conductance curves have been obtained in a range of capacity values. NDR and I(V) resonant 
peak have been observed in the serial configuration due to the different electrostatic coupling. In 
parallel configuration a conduction gap has appeared as some other authors have reported, therefore 
two threshold voltages have been obtained. Each Qd acts like an independent conduction channel and 
their electrostatic interaction is a second order effect. 
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